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Abstract 

The present state of EHD turbulence theory as a possible solar wind 
research tool is surveyed. The theory is statistical, and does not make state- 
ments about individual events. It is unreasonable to expect ever to be able to 
explain individual events with turbulence theory . The ensembles considered 
typically have individual realizations which differ qualitatively, unlike equi- 
librium statistical mechanics. Most of the theory deals with highly symmetric 
situations; most of these symmetries have yet to be tested in the solar wind. 

The applicability of MHD itself to solar wind parameters is highly questionable; 
yet It has no competitors, as a potentially comprehensive dynamical description. 
The purposes of solar wind research require sharper articulation. If they are 
to understand radial turbulent plasma flows from spheres, laboratory experiments 
and numerical solution of equations of motion may be a cheap alternative to 
spacecraft. If "real life" information is demanded, multiple spacecraft with 
variable separation may be necessary to go further. The principal emphasis in 
the theory so far has been on spectral behavior for spatial covariances in wave 
number space. There is no respectable theory of these for hi ghl y anisotropic 
situations. A rather slow development of theory acts as a brake on justifiable 
measurement, at this point. 
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I. INTRODUCTION 


One question solar wind research may ultimately have to answer for 
itself is whether it will resemble geography or dynamic meteorology more closely. 
Geography miiist concern itself with events which are specific and in a very real 
sense accidental. Its use of analytical mathematics is limited and may some- 
times be simply decorative. Dynamic meteorology proceeds from the assimiption 
that a largely complete mathematical description can be formd and should be 
pursued, even if the complexity of the differential eqmtions and the incomplete- 
ness of the boundary data guarantee that the program will be a long time coming 
to completion (see, e.g. , Pedlosky, 1979). 

Much of what happens in any individual rainstorm is extraordinarily 
striking, but it cannot usefully be considered in the light of a detailed 
mathematical theory. It will never happen in quite the same way again. Austere 
discipline is required to focus on those aspects of the weather which are at 
least statistically reproducible, and therefore susceptible to a mathematical 
theory . 

Space physics has frequently taken the "event" as its unit of concern. 
Such-and-such a set of fluctuating field signals were seen on such-and-such a 
detector on such-and-such a day. Plausible hypotheses about what might have 
been responsible for the signals are produced, and are buttressed by such mathe- 
matics as lies ready to hand. The kind of botindary and initial data that would 
be necessary to extract shaxp conclusions from the mathematics are invariably 
lacking, and the machineiy for extracting the concliosions often also does not 
exist. A rather subjective opinion is usually necessary at the end as to 
whether or not the "event" has been satisfactorily "explained". This paradigm 
is by now deeply ingrained, and is an \mconscious ingredient in the evaluation 
of many of the papers, say, which one finds in Joxirnal o^ Geophysical Research . 

The subject is at a natural stage to begin to ask what the possibili- 
ties are for making it into a mathematically tighter and more intellectually 
crisp area of endeavor. It is equally natural to inquire into the range of 
available models which have been pursued in comparable and more highly developed 
continuous-media situations such as meteorology, say, or ocenography. 

The purpose of the following material is a consideration of the ade- 
qmcy of the available solar-wind mathematical description to the task of 
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providing a comprehensive dynamical description. If precedents from nearby 
subjects are any guide, there would seem to be only one serioxas contender as 
a model for what such a mathematical description might look like. That is 
classical Navier— Stokes hydrodynamics, which is the basis for such theories as 
there are for the dynamics of the earth's atmospheres and oceans. Some exposure 
to hydrodynamic theory will be assumed — sufficient, at least, to take for 
granted the unquestioned role hydrodynamics plays as theoretical research tool 
in those subjects. The following pages are a survey of the present status of 
magnetohydrodynamic (hereafter: MHD) turbulence theory and its adequacy as an 
off-the-shelf research tool for describing solar wind measurements. The con- 
clusions are not all rosy, and the analysis of the available mathematical 
descriptions and techniques leads to the belief that they should only be applied 
to solar wind data with extreme caution, and perhaps with a sense of humor. 

In order of ascending complexity , the possible dynamical descriptions 
for the solar-wind plasma are: (l) one-fluid magnetohydrodynamics (MHD) ; 

(2) multi— species , charged— fluid hydrodynamics with assumed closures for the 
pressure tensors (equations of state); (3) the Vlasov description in terms of 
particle distribution functions; (k) Vlasov equations modified by adding Fokker- 
Planck collision terms on the right hand sides. Specialized models, such as 
the Chew-Goldberger-Low approximation, which rather arbitrarily drops heat flow 
along magnetic field lines, can be accommodated in various niches in the above 
list. 

If the expected dynamics of the system were linear and non-turbulent , 
at least the first three models could be taken serioiisly as contenders. The 
controlled fiision (CTR) community has done so, gambling on the hope that labora- 
tory experimentalists will be able to produce confined plasmas whose dynamics 
remain linear and at most weakly turbulent. But by anyone's definitions, the 
solar wind's behavior is \anmistakably turbulent and nonlinear. The fluctuating 
magnetic fields, flow velocities, and electric fields are as large as anything 
that can be defined as averages in the local zero-momentinn frame. The time 
history of any comjxDnent of the fields behaves for all practical purposes like 
a random variable. This is the definition of "strong turb\ilence" , if one is 
needed. 

Becaiise it is the only one of the fovtr descriptions that is close 
to being manageable, even numerically, MHD assumes the role of the only serious 


109 



contender for a "strong turbulence" mathematical description. It is far simpler, 
mathematically, than any of the others, and yet the number of strong turbulence 
problems that we can handle with it will be seen to be extremely limited, even 
assuming its correctness. 

The following material is intended as a brief look at nonlinear MHD 
turbulence theory, as it may be considered as a potential solar— wind research 
tool. Section II deals with the applicability of MHD itself to a plasma with 
solar-wind densities and temperatures, stressing the roles of incompressibility, 
collisionality , and the proper analytical form for the essential dissipative 
terms. Section III simimarizes the status of incompressible MHD turbulence 
theory as it has been developed so far, emphasizing the hi^ degrees of 
symmetry required if even the crudest theories are to have extractible conse- 
quences. Section IV suggests some tentative implications of Sections II and 
III for solar -wind research. 

Anticipating the conclusions, one of them is that there is presently 
available at best only an outline of a theoretical framework in which kinds of 
solar wind data that have been collected could be sensibly interpreted. If we 
are serious about wanting to go beyond a largely descriptive understanding of 
the solar wind, a far higher fraction of our effort will have to go into mder- 
standing the basic plasma physics of the mediimi. The analytical and numerical 
tools now in hand are not adequate to the demands being placed on them by the 
sophisticated collection of vast quantities of data, whose quality is far 
higher than any framework available for maJting use of it . 
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II. THE APPLICABILITY OF MHD 

The one-fluid MHD equations, in the simplest form in which they might 
he considered realistic, is 

V-y = 0 ( 1 ) 

P(H + y.Vy) = -Vp + 2^ + pvV^y (2) 

3B o 

^ + y-Vg = B-Vy + nV^B, ( 3 ) 

with a fluid velocity y, and a magnetic field g, a mass density p, a kinematic 
viscosity V, and a magnetic diffusivity n. The solenoidal condition V-B = 0, 
imposed initially, is preserved by Eq. (3). p is the total pressure, magnetic 
plus mechanical, and is obtained from the Poisson equation which results when 
the divergence of Eq. (2) is taken and use is made of Eq. (l) . 

Eqs. (l)-(3) are the simplest closed-form mathematical description 
known for the mechanical motions of a fluid which is both energetic enough and 
electrically-conducting enough to modify the magnetic field imbedded in it. 

Yet the simplicity of Eqs. (l)-(3) is misleading. We are far from being able 
to give analytical solutions except in highly simplified special cases, or in 
the linear limit. For reasons which are by now well known (althoiagh we will 
review them later) , they make demands on computing capability which we cannot 
always expect to meet, even numerically. 

Some of the assumptions which go into the derivation of Eq. (l)-(3) 
are widely known, such as the neglect of the displacement current relative to 
the conduction current, or the assumption that electrostatic forces are capable 
of keeping the electron and ion charge densities approximately equal. Three 
assumptions need to be singled out for mention in connection with the solar 
wind. They are not obviously fulfilled by solar wind parameters, and the 
serious failure of any one of them can leave us with a mathematical description 
which IS even far less tractable than Eqs. (l)-(3). They are: (l) incompressi- 

bility (V.y = 0); (2) scalar dissipation coefficients v and n; and (3) collision- 
dominated inequalities required in the derivation of Eqs. (l)-(3). 
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(1) Incompressibility (V-y = O) 


Incompressibility is an undisputed featvire of normal fluid mechanics 
that is difficult to justify rigorously. It is usually done (Landau and 
Lifshitz, 1959; Batchelor, 1967) by using estimates for the dominant force 
terms in the equation of motion and their effect, through the compressibility, 
on the density p of a moving fluid element. The change in density Ap for a 
fluid element which experiences a change in pressure Ap may be taken to be: 


Ap 


Ap An 
dP/dp 


(4) 


2 _ 

If the medium obeys an equation of state p = p(p). The sound speed is Cg = dP/dp. 
Ap may be estimated by using either the py*Vy term in Eq. (2) or the B*VB/4 tt 
term. (These are expected to dominate the p3y/3t term and the viscous term 
pvV^v in cases which have significant amounts of turbulence.) For V , we will 

use L~^, where L is a characteristic length over which the fields vary. Ap 

2 

may be estimated from the convective term, first, as of order ~pv . In this 
case, the fractional variation in density is small for a typical fluid element 
if 

M ~ Pl!^ = v^/c ^ « 1 (5) 

P PCg s 


as in ordinary hydrodynamics* Then we may eshimahe Ap from the magneiic force 
term as Ap ~ B^/4 it, and instead of (5), we get 


( 6 ) 

p 47TpC 


or that the magnetic pressure shall be small compared to the mechanical pres- 
sure (3 » 1» in conventional plasma physics jargon). If there is a strong 
mean field Bq present which is large compared to the fluctuating B, 

Ap ~ B^B/Utt, and (6) is replaced by 


P 


B" 

4ttpc ^ 
s 


J 

B 


:r- « 1 


(T) 


or that (again in conventional plasma terms) B/B^ « 3. 
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Though the fulfilment of conditions (5) and (7) might give some con- 
fidence in (say) the applicability of incompressible MHD for a tokamak plasma, 
no one of the three inequalities (5)-(7) can be said to characterize the solar 
wind. Yet the solar wind, except for occasional shock transitions, often shows 
surprisingly little density variation. From the point of view of considerations 
presently known, this tendency toward incompressibility is still slightly 
mysterious. 


(2) Dissipation Coefficients 

Derivations from first principles lead to far more elaborate dissipa- 
tive terms than those which appear in Eqs. (.!)-( 3). Only those who have 
actually dragged themselves throu^ a Chapman-Enskog calciaation of magneti zed- 
plasma transport coefficients can probably appreciate the fragility of the 
enterprise, but a widely-accepted derivation due to Braginskii (I 965 ) [see 
also: Book, 1980 ] yields a considerably more involved term for the viscous 

dissipation than that given in Eq. (2). Reverting to component notation, 

pvV V. shoiad be replaced by the ion viscosity term -Z ap,^/3x,, where, in 

1=1 J 

a coordinate system with the z-axis along the magnetic field B , 



xy 

P 


xz 


P 

yz 


p 

zz 


- % (W + W ) - \ (V - W ) - n W 


XX yy 


XX yy' 


*3 xy 


- Ik (w_ + w_) + \ (w„.. - V) + n U 


XX yy 


XX yy 


'3 xy 


P = -n,W + '3 (W - W ) 
yx 1 xy ' XX yy' 


"“zx = -^2^xz - 


yz 


^zy = ^^xz 


-n 


w . 

O ZZ 


( 8 ) 


113 


The coefficients in Eq. (8) axe 


n =0.96 

o B 1 1 


Tin = 


_ 3 




Cl 1 




_ 1 _ 1 

2 0) . 2% 

Cl 


( 9 ) 


3kgT^/2 is an ion thermal energy, and = eB/m^c is the proton gyro frequency, 
is an ion collision time and is given by 


^ k/n nXe^ 


( 10 ) 


where m^ is the ion (proton) mass, n is the proton number density, e is the 
proton charge, and X is the Coulomb logarithm, typically 10 to 20. The rate of 
strain tensor W,, is 



3v . 3Vj^ 


2 

3 



V.y 


( 11 ) 


The viscosity coefficients f|^, differ mainly hy the 

numbers of powers of they contain in the denominators. The largest term 

for solar -wind plasmas is A formal estimate of the n^~containing terms 

at lAU, using measured values for the length scales and typical fluctuating 

velocities 5 leads to the conclusion that the ri^ terms are by orders of magnitude 

the largest terms in the equation of motion I There is nothing else in the 

equations of motion that they could be equated to unless the coefficients 

multiplying the r\ were themselves small. These coefficients turn out to be 
o 

linear combinations of V*v and V_l*v_l (here, the subscript "x” applied to a 
vector means the components perpendic\iLar to B) . Only if V*v:^0 and - 0 

does it appear that the viscous terms can do anything but completely overpower 
every other term in the equation. This may be a more convincing argument for 
incompressibility than any that can be given in the conventional way, as in 
the previous subsection. It does, however, leave an additional constraint, 
incompressibility in the plane perpendicular to B, which is not built into 
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Eqs. (l)-( 3 ). The constraint Vj_.Vj ^=0 does come up in the Strauss ( 19 T 6 ; 
Montgomery, 1982) equations of "reduced" MHD, which are appropriate to the case 
of strong externally-imposed dc magnetic field (3 « 1), but its content with- 
out the presence of such an externally-imposed B-field is far from clear. 

Finally, and perhaps most annoyingly, even if the divergences of v 
Yj. small, that does not mean that the terms containing are negligi- 
ble. The jungle of terms involved in Eqs. (8) and ( 9 ) does not lead, by any 
known asymptotic expansion, to a simple diffusion-like viscous term such as the 
last term of Eq. ( 2 ), at the time of this writing. It is possible that we 
will remain in the unpleasant position of settling for the relatively tractible 
pvV y term as a crude model of short-wavelength dissipation, knowing full well 
that it is not an accurate representation. 

( 3 ) Collisionality 

Such expressions as Eqs. ( 8 ) and ( 9 ) are the output of lengthy, 
tedious Chapman-Enskog calculations which begin with a transport (e.g., 
Braginskii, I965) equation with a Fokker-Planck collision term, and iterate 
about a local Maxwell distribution. The expansion parameter, assumed small, 
is the ratio of the mean collision time (x. for ions, x for electrons) to 
the time scale T over which the macroscopic field variables vary, or equi- 
valently, the ratio of mean free paths to macroscopic length scales. In 
the solar wind, these ratios, rather than being « 1 , are »1 if standard esti- 
mates are used for mean free paths and collision times. From one perspective, 
it is astonishing that MHD has any relevance to solar wind phenomena. It has 
been suggested, not unconvincingly, that the Fokker-Planck collision terms 
which are used to compute expressions such as Eqs. (8) and (9), are improper 
because of the observed high level of turbulence in the solar-wind magnetic 
field. Free-flight straight-line trajectories are used in evaluating collision 
integrals and are cut off at a Debye length, and these may be less than appro- 
priate for a particle following a tangled field line. But these are no more 
than suggestions at this point, and what their implied modification of Eqs. (8) 
and (9) might be has not been suggested. 

In summary, there are three respects at least in which the validity 
of incompressible MHD with scalar dissipation coefficients might legitimately 
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be doubted for solar-wind parameters. Yet is is the only contender among 
mathematical descriptions which have so far proved tractable enough to lead 
to any comprehensive theory of turbulent situations. Even then, we shall see 
in the following section that further severe restrictions are necessary in 
order to have concrete results emerge. 
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III. MACaiETOHYDRODYNAMIC TURBULENCE 

If the previously-enumerated reservations about the validity of 
the MHD description are passed over, it may be noticed that a certain amount 
of relatively clean theory of MED turbiilence has emerged in the last two decades. 
The theory relies on certain idealizations that render it less than wholly 
applicable to real-life solar wind conditions. Applicable or not, it constitutes 
the only presently-existing framework in which statements about the solar wind 
can be made which axe more than impressionistic or anecdotal. Virtually all 
of is for the uniform-density (p = const.) case, and the incompressibility 
restriction is important. No significant body of strong turbiolence theory 
exists for compressible fluids, even for ordinary neutral gases, and it would 
be unreasonable to expect MED to yield where the simpler compressible system 
has not. 

Use of the term "strong turbulence" in the preceding paragraph is 
intended to differentiate it from "weak turbulence" theory, which is a perspec- 
tive which has shaped most thinking about nonlinear disordered processes in 
plasmas since about 1962. In weak turbulence theory (e.g., Montgomery, 1977), 
the emphasis is on systems whose dynamics may be considered to be the interac- 
tion of oscillatory normal modes, whose oscillation period is short compared 
to the characteristic time of transfer of excitations from one normal mode to 
another . Our reasons for discoimting the value of weak turbulence theory in 
discussing the solar wind will become apparent when we write Eqs. (l)-(3) in 
appropriate dimensionless imits. 

We first observe that there are at least three physically distinct 
time scales represented in the dynamics described by Bqs. (l)-(3). If we call 
a typical rms flow speed U^ (in a coordinate system moving with the local -mpan 
velocity of the solar wind) , a typical rms magnetic field strength , a 
typical suitably defined mean magnetic field <B>, and a typical length scale 
over which the fields vary 1/k, then these three time scales may be defined 
as follows. There are: (l) the "eddy turnover time" (kUQ)~^ associated with 

the fluid motions (in the solar wind, often B ~U ] ; (2) the "Alfven transit 
time" (k|<B>|/Airp) ; and (3) two dissipative time scales (k^v)”^ and (k^n)”^ 

which may be the same or different, depending upon the magnetic Prandtl niunber 
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v/ri. The situation becomes more complex when we realize that there is not 
one length scale -^l/k, but a whole spectrum of scales, present at any instant, 
and the and may be defined locally in the wavenumber k as well. In the 

short wavelength range (large k) , the dissipative effects may be dominant, 
while at large scales (small k) they may be negligible. There is no sharp 
dividing line where one passes from one regime to another. 

Weak turbulence theory assigns orders of magnitude to its time scales 
of its excitations once and for all, and makes no provision for these to change. 
Its limitations are apparent in any situation in which there are fluxes of 
excitations in k space which move from one regime to another. 

The point is that it is unacceptable to neglect any of the terms 
in Eqs. (l)~(3). It is important to resist the temptation to try to treat a 
limited range of k in dynamical isolation from the rest, making approximations 
there that do not apply elsewhere in k, because of some inequalities which 
obtain locally. Eqs. (l)-(3) are a package, no part of which can be ignored 
without peril. It might be argued, as in Sec. II, that more terms are needed 
in Eqs. (l)-(3) to do justice to the dynamics of the solar wind; if so, then 
the effect is to complicate an already almost prohibitively difficult problem. 

It cannot be argued that terms can be dropped because they may be "small" in 
certain ranges of k. 

For the solar wind, and the Alfv^n speed = B^/ v^Tip are compar- 
able in the zero-momentum frame. ’The coefficients V and n are uncertain for 

reasons already given, and may not even be well-defined. If the Spitzer formu- 

o 

la for the conductivity a is adopted, r] = c /hi\0 , If we use the Braginskii 
to estimate the viscosity, then V = get, in cgs mits, 

a 2.5 X 10^^ sec""^ 

2 

V ~ rij^/p ~ 3 X 10 cm /sec, 

— B 5 

at a number density of n ~ 10 cm and a temperature of lO'^^K. Both U and C. 

6 o A 

are typically 2 or 2.5 x 10° cm/sec, and the most typical length scales have 

been measured to be L ~ 10^^ cm [e.g., Matthaeus and Goldstein, 1982]. 

We rewrite all velocities in units of U = B / Airp, all lengths in 

o o 

units of L, all times in units of L/U , and all magnetic fields in mits of 
B^. The dimensionless version of Eqs. (l)-(3) becomes 
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V • y = 0 


(12) 


^ Y ’ ^ ~ ~^P ■‘■5 ’ ^ ^ (13) 

and 

+ y . VB = B • Vy + ^ V^B. (lU) 

The dimensionless ntunhers R and S are the Reynolds number and magnetic 
Reynolds number, respectively: R = UqL/v and S = C.L/n = U L/r\ (since U~B / Atp, 

here). For the numbers cited, R ~ 10^^ and S 10^3. These large values put 
us far into the regime of high Reynolds number turbulence, which is the domain 
of applicability of such theory as we have (e.g., Batchelor, 1970; Panchev, 1971). 

The picture of high-Reynolds number fluid turbulence which has served 
as a model for the recent development of magnetofluid t\irbulence theory is due 
to Richardson, G. I. Taylor, and Kolmogoroff, and is elegantly summarized and 
developed in the classic monograph by Batchelor (l970). It does not maJ^e 
reference to specific solutions of the dynamical equations, which are regarded 
as irreproducible random variables. Instead, statements are made about ensemble 
averages, indicated by angiilar brackets < >, which are hoped to be relatively 
smooth and reproducible. Thios a measured component of the magnetic field, 
might be divided up into a "mean" plus a "fluctuation" 6B^: 

or B . = <B . > + 6B . , 

1 1 1 

(15) 

V. = <v. > + (5v. 

Ill 

for the velocity field, and so on. 

What the brackets < > mean experimentally is a tricky question. 

Ideally, they should represent ensemble averages over a very large niunber of 
experiments prepared in the same way, based on measurements made after a fixed, 
elapsed time. Even in the laboratory this is difficult, and in the solar wind 
it is out of the question. What must be done is to conjectiare something like 
an ergodic hypothesis, which makes it possible to equate phase space averages 
(or ensemble averages) and time averages. Because there is in the solar wind 
an inevitable relative velocity between the solar wind plasma and the measuring 
instruments, these time averages are really averages over a space-time trajectory, 
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in the zero-momentiim frame. By the time the various symmetries necessary to 
interpret the data have been invoked, one has assumed a certain fraction of 
the consequences that one would, ideally, have liked for the experiment to 
demonstrate. A shaky consistency is often the most conclusive imaginable 
outcome. 

Very nearly all the results so far on MHD turbulence concern the case 
of homogeneous turbulence , for which the statistical properties of the fields 
B^(x, t) , t) are independent of x. One conventionally works in the 

zero-momentum frame, =0. If the direction of the magnetic field is not 

externally constrained in some way by bo\indary conditions, then = 0. The 

quantities of theoretical interest then are mostly derivable from the covariances 


Ry.(r, t) = <v. (x, t) V (x + r, 

J -L J 

R? . ( r , t ) = <B (x, t) B.(x + r, 
1 J 

■t) = (x, t) B (x + r, 

J J- J 


t)> 

( 16 ) 

t)> 

(IT) 

t)> 

( 18 ) 


which, by the assumption of spatial homogeneity, are independent of x. 

Virtually all serious theoretical attempts in both fluid and magneto- 
fluid turbiiLence so far have centered around such quantities as these covari- 

~R 

ances. Attempts to calculate R. . , R. ., R. . from a closed, deterministic dynami- 

ij ij ij 

cal description have displayed great ingenuity and some results, but nothing 

that is of obvious use for explanation of solar wind phenomena, so far. 

Analytical approaches to data have been concerned with the rotation- 

ally isotropic case. In this case, the tensor description of Eqs, (l6)-(l8) 

contracts drastically. The R. .(x, t), for example, may be Fourier-decomposed 

1 J 

as 


where 


R^jCx, t) = /dk S®j(k, t)e^- ~ 

k k *"k^ 

sf (k, t) = E„(k, (19) 

ij ~ B 


with a single scalar variable E^(k, t) determining the evolution of the covari- 
ance. is the energy spectrum, and is related to the rms fluctuating 
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field variable 6B by 


= f E^(k, t)dk. (20) 

8tt 0 ^ 

Eq. (19) does imply rotational Isotropy, and the presence of a finite mean 
<B.> will not in general permit this. Analytical impediments to a deductive 
theory are best illustrated by illustrating the dynamics in a Fourier decomposi- 
tion of V and B over a large cubical box, assuming periodic boundary conditions: 

v(x, t) = y(k, t)e^~'~ 

B(x, t) = B(k, t)e^~‘~ (21) 

If we make up a large column vector whose i^^ element can be the real or 
imaginairy part of any component of any one of the expansion coefficients B(k, t) 
or v(k, t) , the Fourier decomposed structure of Eqs . (l2)-(lU) can be written 
symbolically as [see, e.g., Orszag 1977, or Kraichnan and Montgomery, 198 O] : 


dX. 

— - = C. X.X, - V.X. 

dt jk ijk j k 11 


( 22 ) 


where the coupling coefficients are known, and of a kinematical nature. 

The dissipation coefficients come from the viscous and resistive (linear) 
terms, and generally increase quadratically with increasing wavenumber magnitude 

The essential problem with any analytical approach to Eq. (22) is 
that the nonlinear (quadratic) terms are much larger, throughout most of k-space 
than the linear dissipative ones. No linearization can be justified. The 
inequality is measured by the Reynolds numbers R and S, so that, partic\ilarly 
for the solar wind, the nonlinearity may be expected to be strong. 

Attempts at ensemble avaraging moments of the X^ lead to an acute 
closxire problem exemplified by (e.g., Orszag, 1977): 


<X.^/2> +V.<X.^> 
dt i 11 


= E 
Jk 




(23) 
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with a corresponding equation for the time derivative of each nth moment in 
terms of the (n + l)st. 

The situation is reminiscent of the BBGKY hierachy derived from the 
Liouville equation, with the difference that no small parameters suggest them- 
selves as bases for perturbation e^ansions. Great ingenuity has been brought 
to bear, particularly by R. H. Kraichnan (1959, 1964, 1975), on the problem of 
closure approximations for the moment hierarchy derived from Eqs . (22), The 
calcialations are lengthy, require (Kraichnan, 1964) extensive numerical analysis, 
and so far have been limited to the isotropic case. Their generalization to 
anisotropic cases poses formidable problems, and has not been done. 

Eq. (23) e:xpresses the growth or decay of the energy in a particular 
Fourier mode as a sum of a large number of contributions from interacting triads 
of modes whose wave numbers sum to zero. Physical intuition is of limited 
utility in assessing the cumulative effect of the large number of these terms 
which contribute to each mode: the expansion in Fourier series (or other 

orthogonal functions) leave behind any simple resolution into forces and responses, 
"frozen-in field lines", or any of the other readily visualizable but often non- 
quantitative conceptualizations in terms of which MHD has often been discussed. 

The C. , or modal interaction coefficients, are smoothly- varying functions of 

1 JK 

wavenumber where they are non-zero. 

The statistical mechanics of the system (22) with all the dissipation 
coefficients set ^ 0 is tractable. In the cases investigated ( Navier-Stokes 
and MHD in two and three dimensions), truncation at a large but finite number of 
expansion coefficients and equations has led to systems which seem to be ergodic . 
Time averages of phase fixnctions are predictable as ensemble averages (canoni- 
cal or micro canonical) based on the constancy of those invariants which are 
still invariant after the truncation. These conclusions have been repeatedly 
verified numerically [Seyler, et al, 1975; Fyfe et al, 1977a, b; Kells and 
Orszag, 1978 ], and they need only to be alluded to here. 

The difficiilty is that the dissipative terms, if non- zero (v^ 5 ^ O) , 
modify the dynamics qualitatively. Even though they may be relatively small 
over a good part of the wavenumber space, they in effect "pull the plug" at 

2 

the high end of wavenumber space. Because they originate from terms like vV v 
2 

and nV B, they become arbitrarily large, when Fourier-represented, at the large 


122 



values of k. The effect of the (conservative) nonlinear terms is basically 
to scramble, in virtually a stochastic way, excitations from one value of k to 
another. Those excitations that find themselves at large values of [kj get 
gobbled up by dissipation. The flow in k space tends to be toward those regions 
which are deficient, relative to the predictions of the non-dissipative equili- 
brium ensembles. The nonlinear scrambling terms continually try to replenish 
the excitations which are being drained away at high k. Raising the Reynolds 
numbers R and S in Eqs , (l3)^(l^) only increases the "dissipation wave number , 
at which the dissipation sets in, but does not make it go away • The prevailing 
opinion is that the integrated dissipation rate for Eqs. (13) s (l^) remains 
finite even as R and S This gives transfer from one part of the wave- 

number spectrum to another a central iX)le in the dynamics that it does not 
have in linear, or nearly linear, systems. 

Very large numbers of Foiirier modes are required to resolve all the 
dynamically important spatial scales, as R, S become large. This provides 
severe limits on numerical attempts to solve Eqs. (13) and (lU) . A pessimistic 
rule of thumb is that one grid point (or finite element, or Fourier coefficient) 
per dimension per unit Reynolds number is required. Thus, a three-dimensional 
(6U)^ simulation (which will not quite fit in core on a CRAY-l) would be required 
to resolve turbulence with a Reynolds number of 6U. This requirement can be 

relaxed somewhat, but not by an order of magnitude a Reynolds number 1000 run 

could probably not ever be resolved on a (64)^ grid, if the Reynolds niomber 
were to be based on the mean length scale in the flow. When one begins to talk 
about Reynolds numbers many orders of magnitude larger , the real limitations 
of foreseeable computers, in dealing with turbulence, become apparent. 

Dimensional analysis, applied to isotropic, homogeneous situations, 
have led to predictions of power laws in wave number space for the energy spectra 
Eg(k), E^(k) in different situations. The predictions differ from fluids to 
magneto fluids , and from two to three dimensions. They are virtually the only 
simple, testable analytic predictions that four decades of turbulence theory 
have been able to come up with. There are ingenious closures of the hierarchy 
of which Eq. (23) is the first member, but they are not simple, and so far 
they all assume higher degrees of symmetry than the solar wind has been shown 
to possess. These dimensional analysis arguments can be grouped under the 
rubric of "cascade theory". 
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High Reynolds mimber turbulence theories and computations have 
been formalized around two general classes of situations; "forced" and "decay" 
situations. These are at best loose approximations to actual physical cases, 
but they are as close to a universal or situation-independent problem as can 
be isolated- Because the magnetofluid equations are dissipative, a turbulent 
field will eventually decay away, and one can seek features of the decay which 
may be insensitive to initial spectra. Steady-state situations require a 
source of excitations, or "forcing", that is balanced against the dissipation 
rate, averaged over time. The nature of the forcing, often regarded as band- 
limited in wavenumber space, is usually not restricted very specifically, and 
is often modelled by a random function. The search in turbulence theory, as 
elsewhere in physics, is for soluble situations from which a universal, repro- 
ducible, and transferable core of general behavior can be extracted. 

Cascades and Inverse Cascades 

Power laws and cascade processes are expected for forced situations, 
not for decaying ones, unless there is reason to believe that the lifetime of 
the long-wavelength components is sufficiently great that the short wavelength 
components cannot distinguish them from a maintained "source". Under circum- 
stances that have been discussed at great length in the published literature, 
the following table (Table l) shows what has been done so far in the way of 
conjecturing and establishing inertial subrange exponents for fluids and magneto 
f lui ds . 

There is insufficient scope within this article to review in detail 
the evidence and arguments for and against inertial -sub range power laws which 

have been accumulated. There is little doubt that the question of exponents 
has come to occupy more of the territory than it deserves, to some extent 
because there are concrete theoretical predictions. The exponents derive not 
from any dynamical arguments but from conjectured similarity variables. Deriv- 
ing them from dynamics has been the most pursued of all subjects in turbulence 
theory, but no wholly satisfactory resolution has been achieved. Even if it 
were achieved, relatively little light would be shed on the dynamics of the 
solar wind. 
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Table 1 


Situation 

NAVIER -STOKES , 
3D 

NAVIER-STOKES , 
2D 

MHD, 

3D 

MHD, 

2D 

Cascaded 

Quantity 

ENERGY 

ENERGY & 
ENSTROPHY 

ENERGY & 
MAGNETIC 
HELICITY 

ENERGY & 
MAGNETIC 
POTENTIAL 

Direction 
of Cascade 
in k space 

ENERGY UP 

ENSTROPHY UP, 
ENERGY DOWN 

ENERGY UP, 

HELICITY 

DOWN 

ENERGY UP, 
MAGN. 

POT. DOWN 

Predicted 
Power Law, 
Energy 
Spectrum 

KOLMOGOROFF- 
OBUKHOV (19I+I) 

k"^, 

KRAICHNM (1967), 
BATCHELOR , 

LEITH 

k- 5/3 or 

FRISCH et 
al (1975)* 

k-5/3. 

k-"''3 

FYFE et al 
(1977b) 

Experimental 

Verification 

Attempted 

GRANT, 

STEWART, & 
MOLLIET (1962) 

NO 

MATTHAEUS & 
GOLDSTEIN 

(1982) 

NO 

Computation- 
al Verifica- 
tion 

Attempted 

NO 

FYFE et al, 

(1977) 

[insufficient 
spatial resolu- 
tion] 

NO* 

FYFE et al 
(19TT) 

[insufficient 
spatial resolu- 
tion] 


* See also Meneguzzi et al (1981). 


Table 1 . Cascades, Inverse Cascades, and Power Law Predictions. Original 
references are cited in bibliography . 
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The strongest single limitation which the present cascade and inverse 
cascade theory may have when applied to the solar wind concerns the assimiption 
of isotropy, which underlies all of the predictions listed in Table 1, and all 
of the dimensional analysis arguments formulated in the Kolmogoroff style since 
19^1- We are virtually certain that the solar wind is not isotropic, and the 
weaker assumption of axisymmetry may be regarded as open to serious legitimate 
doubt. It is naive to regard the removal of the isotropic restrictions on 
cascade power-law predictions as only a technical point which is sure to be 
overcome soon; its status is at present very dark, and no resolution is in sight 

Selective Decays 

An even more tentative class of generalizations, not without implica- 
tions for the solar wind, are those processes called selective decays , in which 
all the fields decay as in the initial value problem, but some of the global, 
non-dissipative invariants may decay less rapidly than others. Qualitatively, 
there are two possible reasons for this. First, the dissipation is effective 
only at the shorter wavelengths, and quantities transferred to long wavelengths 
may simply stay out of reach of the dissipation. Second, dissipation integrands 
for some variables may be peaked at higher wavenumbers than for others and to 
be more effective at dissipation for this reason. Arguments and computations 
for these possible "selective decay" processes have been presented by Montgomery 
Turner and Vahala (1978), by Matthaeus and Montgomery ( 1980 ), and by Riyopoulos, 
Bondeson, and Montgomery (.1982). 

Each such selective decay process, if valid, would imply a tempor- 
ally decreasing magnitude of the ratio of two of the ideal invariants: energy 

to magnetic helicity for 3D MHD, for example (Taylor 197^ made use of such an 
assumption in predicting asymptotic states of decaying toroidal Z pinches). A 
variational problem arises by minimizing this ratio, which often has for its 
solution a relatively simple Euler equation which predicts a quiescent state. 
Needless to say, this is an attractive possibility. If the tendency of highly 
disordered turbulent motions is to decay to some universal non-trivial quies- 
cent state, regardless of the path of the decay, then this is indeed a wonder- 
fully simple ingredient to add to the few pieces of general information we 
have about turbulence. 
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For example, if for 3D MD, the energy-to-helicity ratio were to 

decay toward its minimimi value, this is simply a force-free state, a solution 

of V X B = AB, where A is a Lagrange multiplier, and y = 0 everywhere. For 

2D MHD, the decay of energy to mean square magnetic potential again leads to 

a quiescent state with a mean magnetic field derivable from a vector potential 

B=VxA=VxAe , where 

z’ 

(V^ +X ^)A = 0. 

Some numerical evidence has been presented for both of these kinds of selective 
decays. The results are encouraging but shoiild be held with extreme caution; 
the conclusions are difficult to document and expensive, and very few runs have 
been carried out [Matthaeus and Montgomery I 98 O; Riyopoulos et al 1982 ]. 

A second class of decay hypothesis, not entirely consistent with the 
first, concerns the ratio of the cross-helicity (another ideal invariant) to 
the energy. The ratio 

= 2 -^y*D dx 

c /(Y^+B2)dx 

a constant in the absence of dissipation, has been shown under some circum- 
stances to increase monotonically (Grappin et al, 1982; Matthaeus, Goldstein, 
and Montgomery 1982) with time in the presence of dissipation. This increase 
points to an equipartitioned state, certainly not quiescent, with y = ±B. 

From the point of view of solar wind observations, this is an attractive 
possibility, because many observations, from Belcher and Davis (l97l) on, 
have shown solar wind velocity fields and magnetic fields to be closely aligned 
or anti-aligned. These are sometimes referred to as ”Alfvenic fluctuations". 

The paradox of MHD turbulence’s tending apparently both to states 
in which is maximal and helicity to energy is also maximal is an example 
of the wide-open character of research into MHD turbulence. There is compelling 
evidence for both conjectures, but both cannot be simultaneously true. If 
either is true, it may well determine the asymptotic state toward which solar 
wind turbulence is trying to decay. 
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IV. SUMMARY 


MHD turbulence theory provides the most nearly adequate framework 
in which to discuss the physics of solar wind turbialence. The collected data, 
however, are far superior both to the available jxzstification of the MHD des- 
cription and to its systematic development for turbiilent fields which lack 
high degrees of symmetry such as rotational isotropy. Ikiderstanding the physics 
of the solar wind at the present time is probably more limited by the unanswered 
questions in turbulence theory than by any scarcity of measurements. Expanded 
experimental programs to probe solar wind txirbxHence, such as that advocated by 
the 1980 Plasma Turbulence Explorer Panel (Montgomery et al, I980) would require 
a considerably broader attack than has so far been mounted on the basic plasma 
physics of the turbulent medium. 
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